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Abstract 

This paper is a sequel to the author's paper entitled "On Dark Matter, Spiral Galaxies, and the 
Axioms of General Relativity" which explored a geometrically natural axiomatic definition for dark 
matter modeled by a scalar field satisfying the Einstein-Klein-Gordon wave equations which, after 
much calculation, was shown to be consistent with the observed spiral and barred spiral patterns 
in disk galaxies, as seen in Figures 3, 4, 5, 6. We give an update on where things stand on this 
"wave dark matter" model of dark matter (aka scalar field dark matter and boson stars), an interesting 
alternative to the WIMP model of dark matter, and discuss how it has the potential to help explain the 
long-observed interleaved shell patterns, also known as ripples, in the images of elhptical galaxies. 

In section 1, we begin with a discussion of dark matter and how the wave dark matter model com- 
pares with observations related to dark matter, particularly on the galactic scale. In section 2, we show 
explicitly how wave dark matter shells may occur in the wave dark matter model via approximate solu- 
tions to the Einstein-Klein-Gordon equations. How much these wave dark matter shells (as in Figure 2) 
might contribute to visible shells in elhptical galaxies (as in Figure 1) is an important open question. 

1 Introduction 

What is dark matter? No one really knows exactly, in part because dark matter does not interact signif- 
icantly with light, making it invisible. Hence, everything that is known about dark matter is due to its 
gravity, which is significant since it composes 23% of the mass of the universe [14]. "Regular" matter 
which makes up the periodic table, typically referred to as "baryonic matter," composes only about 5% 
of the mass of the universe [14]. Dark energy, not to be confused with dark matter, is predicted by and 
is one of the great successes of general relativity. It is a very small but positive energy density with neg- 
ative pressure spread evenly throughout the universe in a manner that explains the observed accelerating 
expansion of the universe [34, 33]. Dark energy composes about 73% of the mass of the universe [14]. 
However, since dark energy is evenly spread throughout the universe, it is an insignificant component of 
the mass of galaxies. 

Dark matter, on the other hand, is most of the mass of galaxies [5, 6]. These dark matter "halos" 
in galaxies are roughly spherically symmetric even when the visible baryonic matter is mostly in a disk 
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Figure 1: Shells in Elliptical Galaxies: NGC 474 (left) and NGC 4382 (right) are two examples of 
elliptical galaxies with shells that are prominent. Left Image Credit and Copyright: P.-A. Due (CEA, 
CFHT), Atlas 3D Collaboration. 

[5, 6]. Hence, while spiral galaxies are disk-like in terms of the light they emit, by mass they are roughly 
spherically symmetric. Elliptical galaxies are also dominated by these dark matter halos. For a summary 
of what is known about dark matter, see [35, 4, 15, 39, 30]. 

It may be reasonable to look at detailed images of galaxies to understand dark matter since these are 
dark matter dominated systems. The two main types of galaxies are disk galaxies, most of which are 
spiral galaxies, and elliptical galaxies. The common occurrence of spiral galaxies could be a major hint 
as to the nature of dark matter. The author explored this idea in detail in [7] and produced simulated 
images of disk galaxies with spiral and barred spiral patterns very similar to photos of actual galaxies, as 
can be seen in Figures 3, 4, 5, and 6. 

1.1 Shells in Elliptical Galaxies 

As spectacular as spiral patterns in disk galaxies are, there is also a fascinating phenomenon in the 
images of elliptical galaxies - the common occurrence of interleaved shells, also known as ripples, in 
their images. While shells in elliptical galaxies are usually hard to see with the unaided eye, computer 
enhancement of photos of elliptical galaxies clearly show partial circular steps in the luminosity of the 
elliptical galaxies typically on the order of 3% to 5% of the galaxy's surface brightness profile [5]. These 
shells have a global structure within the galaxy in that the centers of the partial circular shells lie at the 
galactic center. In addition, the radii of the shells on the two opposite sides of the galaxies are typically 
interleaved, in that the radii of the shells in increasing order alternate between two sides of the galaxy 
as explained on page 202 of [5]. These shells, present in 10% to 20% of all elliptical galaxies [5], or 
perhaps 30% to 50% "depending on how closely one looks" [6], have been known to exist since 1980. 
See Figure 1 for two prominent examples of shells in elliptical galaxies. Quoting from [5], page 203: 
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Figure 2: Shells in Elliptical Galaxies: Simulated images of wave dark matter density projected into 
the plane clearly show interleaved shells which result from first degree spherical harmonic perturbations 
from spherically symmetric scalar field solutions to the Klein-Gordon equation. 

In the classical dynamical model of an elliptical, phase space is populated very smoothly. 
Therefore, the existence of ripples directly challenges the classical picture of ellipticals. One 
likely possibility is that ellipticals acquire ripples late in life as a result of accreting material 
from a system within which there are relatively large gradients in phase-space density. Sys- 
tems with large density gradients in phase-space include disk galaxies and dwarf galaxies: 
in a thin disk, the phase-space density of stars peaks strongly around the locations of circular 
orbits, while in a dwarf galaxy all stars move at approximately the systemic velocity, so that 
there is only a small spread in velocity space. 

Numerical simulations suggest that ripples can indeed form when material is accreted 
from either a disk galaxy or a dwarf system - see Barnes and Hernquist (1992) for a review. 
Moreover, simulations have successfully reproduced the interleaved property of ripples de- 
scribed above. Despite these successes, significant uncertainties still surround the ripple 
phenomenon because the available simulations have important limitations, and it is not clear 
how probably their initial conditions are. 

Implicit in most if not all studies of shells and ripples in elliptical galaxies to date is that dark matter 
is correctly modeled as a WIMP (Weakly Interacting Massive Particle). In this paper as in the previous 
paper [7], we study a geometrically natural axiomatic definition for dark matter, called "wave dark mat- 
ter," as an interesting alternative model of dark matter. Furthermore, in section 2 of this paper we show 
how this geometric model for dark matter may form interleaved shells in its galactic density profile, as 
shown in Figure 2. 

We suggest that this global interleaved shell structure in the wave dark matter exemplified in Figure 
2, while invisible itself, could conceivably play a role in the formation of visible interleaved shells in 
some elliptical galaxies, perhaps as exemplified in Figure 1, through gravity, friction, dynamical fric- 
tion, and the processes mentioned in the above quote. This conjecture is based primarily on the striking 
qualitative similarity between the images in Figures 1 and 2, namely that both sets of images have inter- 
leaved shells. The basic idea that regions of increased dark matter density could gravitationally attract 
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Figure 3: NGC 1300 on the left, wave dark matter simulation [7] on the right (showing gas, dust, and 
stars, not dark matter). Left photo credit: Hillary Mathis/NOAO/AURA/NSF. Date: December 24, 2000. 
Telescope: Kitt Peak National Observatory's 2.1-meter telescope. Image created from fifteen images 
taken in the BVR pass-bands. 

more regular visible matter and hence contribute to visible shells is plausible but is also a very subtle 
question to study. Resolving this conjecture with a high degree of confidence may require sophisticated 
computer simulations based on expert modelings of many different astrophysical processes, making this 
a challenging conjecture to study carefully. 

The primary goal of this paper is to show how interleaved shells may occur in approximate wave 
dark matter solutions to the Einstein-Klein-Gordon equations, done in section 2. How much these wave 
dark matter shells (as in Figure 2) might contribute to visible shells in elliptical galaxies (as in Figure 1) 
is an important open question. 

1.2 Geometric Motivation for Wave Dark Matter 

The only reason dark matter is known to exist is because of its gravity, which affects the orbits of stars 
and causes gravitational lensing of light. But according to general relativity, gravity is just an effect of the 
curvature of spacetime. Hence, at this point, dark matter could just as legitimately be called "unexplained 
curvature" on the scale of galaxies and higher. 

To this end, in [7] we reexamined the question of what the fundamental assumptions of general 
relativity are exactly. In doing so, we defined two axioms, called Axiom and Axiom 1 in [7]. At the 
highest level, one could say that 

• Axiom = "Geometry describes the fundamental objects of the universe." 

• Axiom 1 = "Analysis describes the fundamental laws of the universe." 
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Figure 4: NGC 43 14 on the left, wave dark matter simulation [7] on the right (showing gas, dust, and 
stars, not dark matter). Left photo credit: G. Fritz Benedict, Andrew Howell, Inger Jorgensen, David 
Chapell (University of Texas), Jeffery Kenney (Yale University), and Beverly J. Smith (CASA, University 
of Colorado), and NASA. Date: February 1996. Telescope: 30 inch telescope Prime Focus Camera, 
McDonald Observatory. 

Others may have their own ideas about how to take these two provocative statements and turn them 
into precise axioms which result in general relativity and generahzations of general relativity. The author 
played this game and came up with his own axioms in [7] which result in general relativity with a 
cosmological constant (dark energy) and a scalar matter field (dark matter perhaps?) described by the 
Klein-Gordon equation. Regular baryonic matter is not part of the core theory resulting from the axioms 
but may be added in manually. 

It turns out that a scalar field model of dark matter is not ruled out at this time as a viable theory 
of dark matter. In fact, astrophysicists, beginning with a quantum theory motivation, have already been 
studying these equations as a possible model for dark matter [38, 16, 18, 19, 20, 2, 1, 36, 37, 13]. 

This "scalar field dark matter" model also goes by the name "boson stars." In this paper, we will refer 
to this model of dark matter by a new name: "wave dark matter." We prefer "wave dark matter" because 
the scalar field representing the dark matter satisfies the Klein-Gordon equation, a wave equation on the 
spacetime. Hence, the mental picture one should have for the qualitative characteristics for wave dark 
matter are waves on a pond, but with gravity and one more very important difference: Unlike waves 
on a pond which have a single characteristic velocity, group velocities of solutions to the Klein-Gordon 
equation may be anything less than the speed of light, with smaller group velocities for long wavelengths 
and larger group velocities for shorter wavelengths. In addition, these waves may interfere with one 
another both constructively and destructively. Thus, if we are looking for a name for this model of dark 
matter which is descriptive of how it behaves, then "wave dark matter" is a natural choice. 

The following axiom is one way to express the idea that "Geometry describes the fundamental objects 
of the universe." 
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Figure 5: NGC 3310 on the left, wave dark matter simulation [7] on the right (showing gas, dust, and 
stars, not dark matter). Left photo credit: NASA and The Hubble Heritage Team (STScI/AURA). Ac- 
knowledgment: G.R. Meurer and T.M. Heckman (JHU), C. Leitherer, J. Harris and D. Calzetti (STScI), 
and M. Sirianni (JHU). Dates: March 1997 and September 2000. Telescope: Hubble Wide Field Plane- 
tary Camera 2. 

Axiom 0. The fundamental objects of the universe are a smooth spacetime manifold N, which is both 
Hausdorff and second countable, with a smooth metric g of signature ( — and a smooth connection 
V defined on the interior of N. 

A smooth manifold is a Hausdorff space with a complete atlas of smoothly overlapping coordinate 
charts [29]. Hence, we are reminded that coordinate charts are more than convenient places to do calcu- 
lations, but are in fact a necessary part of the definition of a smooth manifold. Given a fixed coordinate 
chart, let {di}, < z < 3, be the tangent vector fields to N corresponding to the standard basis vector 
fields of the coordinate chart. Let gij = g{di,dj) and Tijk = g{VQ^dj,dk), and let 

M = {gij} and C = {Tijk} and M' = {gij^k} and C' = {TijkA 

be the components of the metric and the connection in the coordinate chart and all of the first derivatives 
of these components in the coordinate chart. 

The analyst in all of us knows that given p.d.e.s on a manifold, even tensorial p.d.e.s which do not 
depend on choices of coordinate charts, the existence and regularity theory for those p.d.e.s, which is 
crucial, reduces to understanding them expressed in some choice of coordinate chart. Hence, in the spirit 
of searching for the easiest theories first, one way to guarantee second order p.d.e.s in coordinate charts 
is to require the action to be no worse than first order quadratic in coordinate charts. These analytical 
considerations inspired the following interpretation of "Analysis describes the fundamental laws of the 
universe." 
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Figure 6: NGC 488 on the left, wave dark matter simulation [7] on the right (showing gas, dust, and stars, 
not dark matter). Left photo credit: Johan Knapen and Nik Szymanek. Telescope: Jacobus Kapteyn 
Telescope. B, I, and H-alpha bands. 



Axiom 1. For all coordinate charts ^ : Q C N ^ and open sets U whose closure is compact and in 
the interior ofQ, [g^V) is a critical point of the functional 



with respect to smooth variations of the metric and connection compactly supported in U, for some fixed 
quadratic functional QuadM with coefficients in M. 

Note that we have not specified the action, only the form of the action. As is standard, we define 



for some functions {F""} to be a quadratic expression of the {xa} with coefficients in Y. 

The implications of the above axioms when the connection is assumed to be the Levi-Civita connec- 
tion have been long understood. When the integrand in Axiom 1 is reduced to QuadM{M'), vacuum 
general relativitiy generically results. When the integrand is reduced to QuadM {M' U M), vacuum gen- 
eral relativity with a cosmological constant generically results. Here "generically" means for a generic 
choice of quadratic functional, so that the zero quadratic function, for example, is not included in these 
claims. These two results were effectively proved by Cartan [10], Weyl [41], and Vermeil [40] and pur- 
sued further by Lovelock [22]. The point to keep in mind is that since (5, V) must be a critical point of 
this functional in all coordinate charts, then something geometric, that is, not depending on a particular 
choice of coordinate chart, must result. Hence, if we remove the assumption that the connection is the 
standard Levi-Civita connection, then the above set of axioms seems like an interesting place to start. 





(1) 
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Figure 7: M 87 (left) and NGC 1132 (right) are typical examples of the majority of elliptical galaxies 
without visible shells. Left Image Credit and Copyright: ING Archive and Nik Szymanek. Date: 1995. 
Telescope: Jacobus Kapteyn Telescope. Instrument: JAG CCD Camera. Detector: Tek. Filters B, 
V, and R. Right Image Credit and Copyright: NASA, ESA, and the Hubble Heritage (STScI/AURA)- 
ESA/Hubble Collaboration. Acknowledgment: M. West (ESQ, Chile). 

We note that Einstein and Cartan famously played around with removing the assumption that the 
connection was torsion free, while still assuming metric compatibility. However, as our beginning point, 
we make neither assumption. Also, Einstein and Cartan were not trying to describe dark matter and thus 
had different objectives in mind. 

We refer the reader to the appendices of [7] for a detailed discussion of the implications of Axioms 
and 1 . The resulting Euler-Lagrange equations for certain actions satisfying these axioms are, quite 
remarkably, equivalent to the Einstein- Klein-Gordon equations with a cosmological constant A: 



for a real scalar field /, where / = everywhere corresponds to the Levi-Civita connection [7] and 
there is a formula for the connection of the manifold in terms of f and the Levi-Civita connection. Note 
that / = recovers vacuum general relativity with a cosmological constant. In the above equations, the 
speed of light and the universal gravitational constant have been set to one, G is the Einstein curvature 
tensor, □ is the wave operator (the divergence of the gradient) of the Lorentzian metric g, f is the scalar 
field which we suggest causes the curvature of spacetime attributed to dark matter via the right hand side 
of equation 3, A is the cosmological constant, and T is another (new) fundamental constant of nature. 
Estimates on the value of T are computed in [8] by comparing spherically symmetric solutions of the 
Einstein-Klein-Gordon equations to data collected about dwarf spheroidal galaxies. 
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Figure 8: Simulated images of the same wave dark matter density as in Figure 2 projected into the plane 
but nearly along the axis of the shells hides the interleaved shells which result from first degree spherical 
harmonic perturbations from spherically symmetric scalar field solutions to the Klein-Gordon equation. 

It is more typical for the Klein-Gordon equation to describe a complex scalar field /. If one is will- 
ing to allow the connection coefficients Fjjfc to be complex valued, then complex valued / are natural 
as well. For the Einstein-Klein-Gordon system, a complex scalar field is equivalent to two real scalar 
fields. Complex scalar fields are more convenient because they have static spacetime solutions [32, 31]. 
Analogous "nearly static" real solutions can be achieved as \/2Re(/) whose spacetime metrics are very 
similar except for a high frequency oscillating pressure term which averages out to zero. Preliminary es- 
timates suggest that this oscillating pressure effect might not be significant enough to make it measurable 
on physically relevant time scales. Hence, there may not be a testable difference between the predictions 
of real and complex scalar field dark matter at this time, though this deserves further thought. In either 
case we will call these scalar field models "wave dark matter." 

Much more is known about the predictions of WIMP dark matter than wave dark matter. WIMP 
dark matter has been extremely successful on cosmological scales [11]. On the other hand, [7] proves 
that wave dark matter, unlike WIMP dark matter, is automatically cold, as observed, if the universe is 
assumed to be homogeneous and isotropic. 

There are many basic questions still open about wave dark matter. Some of the most important ques- 
tions center around understanding which wave dark matter solutions are physically relevant. Spherically 
symmetric static solutions are easy to compute and so are a natural place to start, which raises questions 
about which of these solutions are stable. Results in [1] indirectly suggest that regular matter could help 
stabilize wave dark matter solutions. Other related results are found in [17]. Of course, the most phys- 
ically relevant wave dark matter solutions could be much more complicated. Computer simulations of 
solutions to the Einstein-Klein-Gordon equations with regular matter included as well are needed to gain 
more insight on these questions of central importance. 

If wave dark matter is assumed to clump on the scale of galaxies for example, then it would be 
reasonable to guess (and important to verify) that it has very similar predictions as WIMP dark matter on 
the cosmological scale. In [23, 24] it is argued that the scalar field dark matter model, which is effectively 
the same as the wave dark matter model discussed here, should reproduce some of the successes of the 
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standard ACDM model above galactic scales. Hence, comparing and contrasting the predictions of 
WIMP dark matter and wave dark matter on the scale of galaxies might be required to find the most 
obvious differences in their predictions. 

1.3 Dark Matter on the Galactic Scale 

In this subsection we list six conclusions about dark matter based on astronomical observations. We 
have roughly ordered the list by confidence in the statement and have placed parentheses around the 
statements which are less certain. 

Conclusions from Astronomical Observations 

1. Disk galaxies commonly have spiral and barred spiral density wave patterns in their images [5, 6]. 

2. Elliptical galaxies commonly have interleaved shells in their images [5,6]. 

3. Most of the dark matter in the universe is cold (traveling at non-relativistic speeds) [5, 6]. 

4. (There is a rough lower bound on the mass of isolated "blobs" of dark matter [9, 35].) 

5. ((The dark matter density in the cores of galaxies appears to be bounded [35, 21, 2, 12].)) 

6. (((It is possible that there are oscillations in the dark matter density as a function of radius, as in 
Figure 9, for some dwarf spheroidal galaxies.))) 

There are many other known facts about dark matter in galaxies that we have not listed here [35]. 
For example, the rotation curves of spiral galaxies provide a tremendous amount of information about 
the distribution of dark matter in spiral galaxies and suggest the notion of a universal rotation curve [35]. 
Also, the TuUy-Fisher relation [25, 26], that the luminous mass of a spiral galaxy divided by the typical 
velocity of stars in the galaxy to the 4th power is roughly a constant for every spiral galaxy in the known 
universe, is a very remarkable fact deserving of a great deal of further study and consideration. Either of 
these last two observations could be extremely important for revealing the true nature of dark matter. 

For now, though, we focus on the six conclusions listed above because plausible statements can be 
made differentiating wave dark matter and WIMP dark matter with regards to these conclusions based 
on astronomical observations. The first two conclusions are observations and hence are beyond dispute. 
The third conclusion is for all intents and purposes beyond dispute given that dark matter appears to be 
gravitationally trapped in galaxies and clusters of galaxies, which would not be the case if it were moving 
at speeds close to the speed of light. 

The last three conclusions, however, are somewhere between probable and possible. Naturally it is 
frustrating not to have complete clarity on important issues, but given that dark matter is invisible and 
known only by its gravity, we have to take what we can get. 

The smallest known "blobs" of dark matter are dwarf speroidal galaxies. In addition, dwarf spheroidal 
galaxies are almost entirely dark matter [35], which is one reason they may be one of the best sources of 
information about dark matter. Another curious fact is that there appears to be a rough lower bound on 
the masses of dwarf spheroidal galaxies, below which "blobs" of dark matter have not been observed to 
exist [9, 35]. By rough, we include the possibility that the frequency of such blobs of dark matter drops 
off quickly as a function of total mass at some point. On the other hand, we have put conclusion 4 in 
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Figure 9: Density as a function of radius for six dwarf spheroidal galaxies based on the distribution of 
the velocities of the stars in each galaxy [42]. Note that three densities have oscillations in their density 
profiles (arbitrarily dashed in gray starting at the first local minimum). Given the uncertainties in the 
data and the assumptions of the model, it is unclear whether or not these oscillations in the dark matter 
density actually exist. 

parentheses since dark matter is, after all, invisible. Hence, maybe smaller blobs of dark matter exist, but 
we just can not see them. This issue is commonly referred to as "The Missing Satellites Problem" [9]. 

When WIMP dark matter is modeled on computers under a certain range of initial conditions, it is 
found to settle down to a density distribution which is well approximated by the spherically symmetric 
Navarro-Frenk- White profile given by 



for some given values of pQ and the scale radius Rs [27, 28, 12]. In particular, the density at the origin 
(assumed to be the center of mass of the system) diverges to infinity in this density profile like poRg/r. 
Given this, astronomers have looked for evidence for the dark matter density in galaxies to be unbounded 
at their centers. However, so far no such density spikes in dark matter have been found, which is the 
content of conclusion 5, sometimes referred to as "The Core-Cusp Problem" [12]. 

The weaknesses of conclusions 4 and 5 is that they concern "not observing" something as opposed 
to directly observing something. We have ranked conclusion 4 higher than conclusion 5 because masses 
of galaxies are found in a wide continuum right down to the lower limit at which point there is close to 
nothing. Hence, the actual observation of galaxies with masses in the entire range above the cutoff is 
very relevant. Nevertheless, concerning Conclusion 5, astronomers have had better success fitting what 
they call "cored" density profiles (Navarro-Frenk- White with the spike removed, like coring an apple) to 
dark matter densities deduced in galaxies from observations. 




(5) 
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Finally, conclusion 6 is the least certain conclusion of the bunch and falls into the "possible" category. 
More observations and analyses are required to determine whether observation 6 is really true, so the 
reader is duly cautioned. Nevertheless, conclusion 6, if true, could have far reaching implications about 
the nature of dark matter, so this author decided it deserved mentioning. 

Conclusion 6 is based entirely on the work in [42], though those authors do not claim this conclusion 
in their paper. In [42], the authors attempt to compute the spherically symmetric density as a function of 
radius of six dwarf spheroidal galaxies based on the observed velocities of the stars in each galaxy. This 
type of inverse calculation has a list of assumptions built into it. In addition, there are error bars in such 
calculations coming from uncertainties in the data which were not clear to this author from reading the 
paper. Nevertheless, the result of the model used in [42] is that half of the six dwarf spheroidal galaxies 
had oscillations in their densities as a function of radius in order to best fit the data, as seen in Figure 9. 
Given the uncertainties in the data and the assumptions of the model, it is unclear to this author whether 
or not these oscillations in the dark matter density in these three dwarf spheroidal galaxies can actually be 
concluded to exist. However, even if the error bars are such that these oscillations may not be concluded 
to exist, one may conclude that oscillations in dark matter density as a function of radius are at least 
compatible with these particular observations. 

In fact, in [42], the authors suggest these dark matter density oscillations are not real, though there is 
very little discussion. Given the pervasiveness of the WIMP model of dark matter and the fact that these 
oscillations are highly unexpected from a WIMP dark matter point of view, this author was concerned 
that the possibility of dark matter density oscillations in dwarf spheroidal galaxies may have been dis- 
carded prematurely. Hence, we have included conclusion 6, perhaps better named "possibility 6," here 
to stimulate more discussion of this provocative possibility. 

1.4 Wave Dark IVIatter on the Galactic Scale 

So how does wave dark matter compare to WIMP dark matter in terms of predictions each makes on 
the galactic scale? Our discussion here, while relevant, does not come close to setthng this complex 
question. 

In this subsection, we discuss the ways in which the six "Conclusions from Astronomical Obser- 
vations" about dark matter in galaxies from the previous subsection favor the possibility of wave dark 
matter as a viable dark matter model. We believe that the discussion presented here justifies wave dark 
matter as a legitimate dark matter candidate well deserving of further study. 

1. Disk galaxies commonly have spiral and barred spiral density wave patterns in their images 
[5, 6]. 

2. Elliptical galaxies commonly have interleaved shells in their images [5, 6]. 

Even though most of the mass of galaxies is dark matter, in the standard explanations of these two 
phenomena dark matter plays a secondary role. If this is true, wave dark matter and WIMP dark matter 
are both able to fit this picture since both have spherically symmetric solutions. However, looking at 
first and second degree spherical harmonic perturbations from spherical symmetry for wave dark matter 
solutions gives an additional possible mechanism by which to explain each of these phenomena. In these 
wave dark matter mechanisms, the substructure of the dark matter (the perturbations from spherical 
symmetry) is what drives the substructure in the regular matter (spiral patterns and shells). 
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In the standard Lin-Shu density wave theory, spiral galaxies are modeled with spiral density waves 
resulting from global gravitational instabilities in the visible matter. The dark matter is modeled as 
spherically symmetric and any dark matter substructure is assumed to be of secondary importance. By 
contract, in [7] we show how second degree spherical harmonic deviations from spherical symmetry in 
the wave dark matter can lead to rotating ellipsoidal galactic potentials which can produce spiral patterns 
in the gas, dust, and stars, thereby accounting for both spiral and barred spiral patterns in disk galaxies, 
as demonstrated in Figures 3, 4, 5, and 6. 

Similarly, the standard explanations for the observed interleaved shells in elliptical galaxies as sum- 
marized on pages 202-203 of [5] also model dark matter as being roughly spherically symmetric and 
make no reference to dark matter substructure. By contrast, in this paper we show how first degree 
spherical harmonic deviations from spherical symmetry in the wave dark matter leads to interleaved 
shells in the dark matter density, as shown in Figure 2. As we suggested in this paper already, this global 
interleaved shell structure in the dark matter, while invisible itself, could conceivably play an important 
role in the formation of visible interleaved shells of regular matter through gravity, friction, dynamical 
friction, and other processes, though this is a complex open question. 

Note that the perturbations in the wave dark matter solutions that we are suggesting are first degree 
spherical harmonic perturbations to explain shells in elliptical galaxies and second degree spherical har- 
monic perturbations to explain spiral and barred spiral patterns in disk galaxies. Our reasoning for this 
difference goes as follows: Typically, we would expect first degree spherical harmonic perturbations to 
be larger than second degree spherical harmonic perturbations which, in turn, we would expect to be 
larger than third degree spherical harmonic perturbations, etc. This picture is compatible with our mod- 
eling of shells in elliptical galaxies, but raises the question of why we assume second degree spherical 
harmonic perturbations (instead of first degree) typically dominate is spiral galaxies. 

Our reasoning is that first degree spherical harmonic perturbations of wave dark matter are the only 
ones which do not fix the center of mass of the wave dark matter at the origin. Since the overall center 
of mass is fixed at the origin, this means that the regular matter center of mass would have to oscillate as 
well. We speculate that spiral galaxies have a roughly fixed center of mass for the regular matter, perhaps 
due to the important role played by friction of the gas and the dust in the galaxy. Hence, a roughly 
fixed center of mass for the regular matter in a spiral galaxy would not allow for a first degree spherical 
harmonic perturbation to dominate, in which case we would expect a second degree spherical harmonic 
perturbation of the wave dark matter to dominate. 

While wave dark matter substructure offers tantalizing possibilities to explain visibile matter sub- 
structure in the images of both spiral and elliptical galaxies, much more study is warranted. Detailed 
computer simulations which model the relevant complicated astrophysical processes correctly are mostly 
likely required to make substantial further progress. 

3. Most of the dark matter in the universe is cold (travehng at non-relativistic speeds) [5, 6]. 

Unlike WIMP dark matter, wave dark matter must be cold in a homogeneous isotropic universe [7]. 
Since dark matter is observed to be cold, it is reasonable to call this a successful prediction of wave dark 
matter. Had the dark matter in the universe been observed to be hot or warm and homogeneous and 
isotropic, this would have ruled out wave dark matter as a reasonable dark matter model. 

The difference in this regard between wave dark matter and WIMP dark matter has to do with the 
wave nature of wave dark matter and the particle nature of WIMP dark matter. In WIMP dark matter, the 
WIMPs can be taken to be uniformly spread over the universe at every point and going in every direction 
at any fixed velocity one likes. Hence, homogeneous isotropic universes can have WIMP dark matter at 
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any temperature, from zero velocity to close to the speed of light. With wave dark matter, however, the 
stress-energy tensor is not additive in this same way. The wave function describing the wave dark matter 
is additive, but the stress-energy tensor is quadratic in the wave function. As shown in [7], homogeneity 
and isotropy of the universe imply that the wave dark matter is very cold at every point as long as the rate 
of change of the Hubble constant is much smaller in norm than T^, the constant from the Klein-Gordon 
equation. More precisely, in Appendix D of [7], we prove the following theorem. 

Theorem 1. Suppose that the spacetime metric is both homogeneous and isotropic, and hence is the 
Friedmann-Lemattre -Robertson-Walker metric —dt^ + a{t)^ds^, where ds^ is any constant curvature 

metric. Iff{t, x) is a real-valued solution to the Klein-Gordon equation Ogf = T^f with a stress-energy 
tensor which is isotropic, then f is solely a function oft. Furthermore, if we let H{t) = a'{t) /a{t) be the 
Hubble constant (which of course is actually a function oft), and p{t) and P{t) be the energy density 
and pressure of the scalar field at each point, then 



P 



p 1 + e 

where 



(6) 



4T2 

and 

fb 1 fb 



j pit) dt, P = j Pit) dt, H' 



J^H'it)fit)^dt 



(8) 



where a, b are two zeros of f (for example, two consecutive zeros). 



In this theorem, the temperature of the universe is approximated by the quantity ^ . Hence, if the 
rate of change of the Hubble constant H'it) is much smaller in norm than T^, then H' must be as well. 
Hence, e and the temperature expression ^ must be close to zero. This is the precise sense in which 
homogeneity and isotropy of the universe imply that wave dark matter must be cold. 

Of course the universe is not perfectly homogeneous and isotropic, so it is an interesting direction 
to take this theorem and try to study what these types of arguments imply when the universe is only 
assumed to be approximately homogeneous and isotropic. Also, a similar result is true when the scalar 
field is allowed to be complex since one complex scalar field is equivalent to two real scalar fields. 

4. (There is a rough lower bound on the mass of isolated "blobs" of dark matter [35].) 

WIMP dark matter simulations predict clumps of dark matter on all scales, but this is not what is 
observed [9] . Instead, clumps of dark matter have been observed on all scales down to a smallest one, 
the mass of dwarf spheroidal galaxies, and then nothing more [35]. 

It is not definitively clear what wave dark matter predicts in this regard. However, static spherically 
symmetric solutions to the Einstein-Klein-Gordon equations which describe wave dark matter come in 
discrete steps [32] as a ground state, first excited state, second excited state, etc. Each of these "solutions" 
is actually a one parameter family of solutions which allows the total mass of the clump of wave dark 
matter to vary. However, if for some reason there is a natural choice of scale for each of these excited 
states, perhaps given by a preferred time frequency for the wave dark matter, then this would predict a 
lowest mass size for a blob of wave dark matter given by the ground state solution. Understanding wave 
dark matter in this regard is a very good direction to pursue. 
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5. ((The dark matter density in the cores of galaxies appears to be bounded [35, 21, 2, 12].)) 



WIMP dark matter, as discussed in the previous subsection, predicts a spike in the dark matter density 
going to infinity at the centers of galaxies [12, 27, 28]. Wave dark matter, on the other hand, does not. 
Detecting such a dark matter spike would most likely rule out wave dark matter as a reasonable dark 
matter candidate. However, as such spikes have not been observed yet, and as astronomers are leaning 
towards the conclusion that they are not there [35, 21,2], this observation, or lack of one, perhaps favors 
wave dark matter at this time. This issue is important for astronomers to continue to study as detection 
of a spike in the dark matter density or a definitive demonstration that such a spike does not exist has 
important implications for the nature of dark matter. 

6. (((It is possible that there are oscillations in the dark matter density as a function of radius 
for some dwarf spheroidal galaxies.))) 

As already discussed in the previous subsection, it is not clear whether this possibility is actually 
true, so gathering more data and studying the existing data for dwarf spheroidal galaxies is warranted. 
However, qualitatively speaking, oscillations in the dark matter density as a function of radius for some 
dwarf spheroidal galaxies is a very reasonable prediction from the point of view of wave dark matter since 
the static spherically symmetric solutions for wave dark matter all have this feature, except for the ground 
state which has monotonically decreasing density as a function of radius [32]. The radial density function 
for WIMP dark matter, on the other hand, would be expected to be closer to the monotonically decreasing 
Navarro-Frenk- White profile in equation 5. Perhaps our main contribution here is to ask astronomers and 
astrophysicists to continue to try to estimate dark matter density in the smallest galaxies, or perhaps those 
a little bit bigger than the smallest galaxies, to see if oscillations in the dark matter density as a function 
of radius are there or not. It is also a good problem to try to figure out what the precise predictions of the 
two dark matter theories are as well. 

In summary, these six conclusions about dark matter based on observations show that there are many 
questions that remain about the true nature of dark matter. As such, both WIMP dark matter and wave 
dark matter deserve much further study to see how their predictions compare to observations, especially 
on the galactic scale. 

2 Wave Dark Matter in Elliptical Galaxies 

Spacetimes associated with galaxies are well approximated by the Minkowski spacetime. In this paper 
we will focus on solutions to the Klein-Gordon equation on the Minkowski spacetime. Solutions in this 
case are qualitatively similar to solutions to the full Einstein-Klein-Gordon equations in many ways, for 
example in spherical symmetry. One main difference is that there are solutions with finite total mass for 
the Einstein-Klein-Gordon equations which are qualitatively similar to solutions in the Minkowski case, 
but only up to a certain radius representing the point at which Einstein-Klein-Gordon equation solutions 
begin to decay exponentially but Klein-Gordon equation solutions on the Minkowski spacetime do not 
(but instead oscillate forever). As explained in section 4 of [7], qualitatively reasonable approximations 
to solutions to the Einstein- Klein-Gordon equations in equations 3 and 4 can be found by simply cutting 
off, outside a large radius, solutions to the Klein-Gordon equation on the Minkowski spacetime. 
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Figure 10: Plots of fujo,o{r) and r • fui,i{r) defined by equation 11 and essential for understanding the 
wave dark matter density in equation 21. The fact that these two oscillating functions go in and out of 
phase is what causes the interleaved shell patterns in the wave dark matter density. Image created by 
running the madab command shells(500, .1, 100000, 60000, 1, -1, 5, 10), available at the author's "Wave 
Dark Matter Web Page" at http://www.math.duke.edu/~bray/darkmatter/darkmatter.html. 
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Cross Sectional Dark Matter Density in XY Plane 




Figure 11: Cross sectional wave dark matter density in the xy plane. The red axis is the x- 
axis and the green axis is the y-axis. Image created by running the matlab command shells(500, 
.1, 100000, 60000, 1, -1, 5, 10), available at the author's "Wave Dark Matter Web Page" at 
http://www.math.duke.edu/~bray/darkmatter/darkmatter.html. 
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Projected Dark Matter Density in XY Plane 
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Figure 12: Projected wave dark matter density in the xy plane. The red axis is the x-axis and the green 
axis is the y-axis. This is what the wave dark matter would look like if it glowed white, as observed 
from a distant observer looking down from the point of view of the z-axis. Image created by running 
the matlab command shells(500, .1, 100000, 60000, 1, -1, 5, 10), available at the author's "Wave Dark 
Matter Web Page" at http://www.math.duke.edu/~bray/darkmatter/darkmatter.html. 
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Cross Sectional Dark Matter Density 




Figure 13: Cross sectional wave dark matter density in a random plane, where the x-axis is red, the 
y-axis is green, and the 2;-axis is blue. Image created by running the matlab command shells(500, 
.1, 100000, 60000, 1, -1, 5, 10), available at the author's "Wave Dark Matter Web Page" at 
http://www.math.duke.edu/~bray/darkmatter/darkmatter.html. 
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Projected Dark Matter Density 
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Figure 14: Projected wave dark matter density in the same plane as the previous figure, where 
the X-axis is red, the y-axis is green, and the z-axis is blue. This is what the wave dark mat- 
ter would look like if it glowed white. Image created by running the matlab command shells(500, 
.1, 100000, 60000, 1, -1, 5, 10), available at the author's "Wave Dark Matter Web Page" at 
http://www.math.duke.edu/~bray/darkmatter/darkmatter.html. 
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Figure 15: Shells in Elliptical Galaxies: NGC 474 is an example of an elliptical galaxy with 
shells. Image Credit and Copyright: P.-A. Due (CEA, CFHT), Atlas 3D Collaboration. From 
http://apod.nasa.gov/apod/apl 10726.html, the NASA Astronomy Picture of the Day: "The multiple lay- 
ers of emission appear strangely complex and unexpected given the relatively featureless appearance of 
the elliptical galaxy in less deep images. The cause of the shells is currently unknown, but possibly 
tidal tails related to debris left over from absorbing numerous small galaxies in the past billion years. 
Alternatively the shells may be like ripples in a pond, where the ongoing collision with the spiral galaxy 
just above NGC 474 is causing density waves to ripple though the galactic giant. Regardless of the ac- 
tual cause, the above image dramatically highlights the increasing consensus that at least some elliptical 
galaxies have formed in the recent past, and that the outer halos of most large galaxies are not really 
smooth but have complexities induced by frequent interactions with - and accretions of - smaller nearby 
galaxies. The halo of our own Milky Way Galaxy is one example of such unexpected complexity. NGC 
474 spans about 250,000 light years and lies about 100 million light years distant toward the constellation 
of the Fish (Pisces)." 
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Figure 16: Projected wave dark matter density, where the x-axis is red, the y-axis is green, and the z-axis 
is blue. This is what the wave dark matter would look like if it glowed white, from 6 different point of 
view. The xy plane, the xz plane, and the yz plane are shown in the top row. Images created by running 
the matlab command shells(500, .1, 100000, 60000, 1, -1, 5, 10), available at the author's "Wave Dark 
Matter Web Page" at http://www.math.duke.edu/~bray/darkmatter/darkmatter.html. 



Using the standard spherical coordinate conventions that 

X = r sin 6 cos (j) 
y = r sin 6 sin (p 
z = r cos 9 

the Klein-Gordon equation on the Minkowski spacetime in spherical coordinates is 

+ ^ + 2 + l i^^ + cote — + ^——]\f-r^f (9) 

\ dt"^ dr^ r dr ^89"^ 89 sin^ 9 d<f>'^ J } 

Solutions can then be expanded in terms of spherical harmonics to get solutions which are linear combi- 
nations of solutions of the form 

/ = Acos{cot) ■ Yn{9,cl)) ■ r" • /^,„(r) (10) 

where 

r 

and Yni9, (j)) is an nth degree spherical harmonic. Note that we require „(0) = but have not specified 
an overall normalization. Naturally, to get a complete basis of solutions we also need to include solutions 
like the one above but where cos(tjt) is replaced by sin(cjt). We will study real solutions in this paper, 
but complex solutions are quite analogous. 
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The form of the Minkowski spacetime solutions that we will study here is 

f = Ao cos{uiot)fojo,o{r) + Ai cos{uiit -<p)sme -r- /wi,i(r) (12) 

where fujo,o{r) and foji,i{r) satisfy equation 11. We note these solutions fall into the form of equation 
10 since both cos (j) sin 6 and sin (/) sin 6 are first degree spherical harmonics, namely the Cartesian coor- 
dinates X and y restricted to the unit sphere in R^. Hence, the above Minkowski spacetime solution is 
the sum of a spherically symmetric solution (degree n = 0) and two degree one solutions. 

The wave dark matter solution represented by equation 12 can be seen to be slowly rotating in the xy 
plane by making the substitution 

a = 4> — {ui — loq) t (13) 
into our expression for / in equations 12 to get 

f = Ao cos{ujot)fojo,o(r) + A-y cos(a;ot - a) sin6l • r • /a,i,i(r). (14) 

Note that to the extent that a. stays fixed in time, then the above solution does not rotate and gives a fixed 
interference pattern since both terms are oscillating in time with the same frequency wq. Hence, we see 
that we get a dark matter interference pattern which is rotating according to the formula 

^0 = (^^1 - wo)i (15) 

with period 

Tdm = • (16) 

U\ - Wo 

Note that a different combination of first degree harmonic terms could produce different dynamics, not 
just rotations. This raises natural questions about observed shells in ellipticals: do they ever rotate 
and, more generally, what are their dynamics? These are potentially very hard questions because the 
time scales involved could perhaps be on the order of pattern periods in spiral galaxies, which can be 
hundreds of millions of years. 

Our next goal is to approximate the wave dark matter density \iwdm due to this scalar field dark 
matter solution and then to expand it in terms of spherical harmonics. From equation 14, we have that 

/ = cos(a;ot) [Ao^o,o('') + ^iCOsasin6' • r • /a;i,i(r)] 



+ sin(a;oi) [^i sin a sin • r • /a;i,i(r)] . (17) 
Hence, by equation 3, the wave dark matter density is 

mDM = ^G{dt,dt) (18) 

OTT 

~ ^ + (-/*' + IV./P)) + (19) 

~ (y) ^^^^ 



where we have assumed a long wavelength solution in the approximation so that spatial derivatives of / 
are much smaller than the time derivatives. We will also assume that ~ T wi, which allows us to 
think of a as being approximately fixed in time. Hence, 
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fJ-WDM ~ [Aof^Qfl{r) + AicosasinO ■ r ■ f^j^^i{r)f + [Aisinasine ■ r ■ f^-^^i{r)f (21) 
= ^o/a;o,o(r-)' + Ajsin^ 9 ■ ■ U,i{rf + 2^o^i cos a sin ^ • r • /u,o,o(r)/u,i,i(r)(22) 

Note that equation 21 only depends on Aq, Ai, and the two functions fuiofl{r) and fuii,i{r) defined 
by equation 11. If we define k'^^ = — T^, then these two function only depend on /cq and ki, which 
are roughly the spatial frequencies of fijjQ,o{r) and foji,i{r) seen in Figure 10. To be physically relevant 
(more specifically to have group velocities for the wave dark matter much less than the speed of Ught), ojq 
and cji must be close to T, resulting in ko and ki being small and wavelengths for fu)Qfi{r) and f^x^iif) 
much greater than 1/T. This observation, combined with precise measurements of the radii of shells in 
elliptical galaxies thought to be due to wave dark matter, might be able to be used to find lower bounds, 
though not very sharp, on the value of T. This is a good project for another paper. 

2.1 Computer Simulations of the Interleaved Shells in the Dark Matter Density 

The wave dark matter density defined in equation 21 has interleaved shells. The easiest way to visualize 
these shells is to simply compute examples and look at the pictures generated by Matlab code, as the 
author has done. The Matlab algorithm "shells.m" is located on the author's "Wave Dark Matter Web 
Page" at http://www.math.duke.edu/~bray/darkmatter/darkmatter.html. All of the computer generated 
interleaved shells in this paper were created by executing the command Une 

shells(500, .1, 100000, 60000, 1, -1, 5, 10); 

In the above command line, 100000 is the maximum radius to which the dark matter density was com- 
puted and 0.1 was the step size used to solve the o. d. e. in equation 11. 60000 is the radius of the view 
displayed in each generated image. 500 is the number of pixels in each direction in each image. 1 and 
—1 are the values for Aq and Ai. When \Aq\ = |^i|, the magnitude of the interleaved shells relative 
to the overall density is maximized by convention of the overall normalizations of fu)(tfl{r) and /a;i.i(?")- 
The program then lets the user specify the number of major shells desired, which in this case is 5 (on 
each side to a radius of 100000, or 3 to a radius of 60000). Finally, 10 is the number of minor shells 
inside each major shell, as can be seen in Figure 11. The matlab code uses the requested number of 
major shells and the requested number of minor shells within each major shell to compute the spatial 
frequencies k^ and ki described above. In the examples in this paper, k^ 0.0033 for a wavelength of 
about 1900 and k\ k, 0.0030 for a wavelength of about 2100. One could take these numbers to be in the 
units of light-years. Since an overall scaling is allowed, however, the individual values are not relevant; 
only the relative values are important. 

Figure 1 1 shows the wave dark matter density in the xy plane. The red axis is the x-axis, the green 
axis is the y-axis, and in other pictures the blue axis is the z- axis. Note that, ironically, we are coloring 
the dark matter white in every computer generated image of interleaved shells in the dark matter density. 
The next image. Figure 12, shows what the dark matter would look like from a distant observer looking 
down on the xy plane from the point of view of the z axis if the dark matter glowed white everywhere. 
These "projected dark matter densities" are created by taking the entire three dimensional wave dark 
matter density and projecting it additively into the chosen plane. Figures 13 and 14 are the same as 
Figures 11 and 12 except that the point of view has been changed. Note that the blue z-axis is now 
visible. Also note that the image in Figure 14 is the same as the second image in Figure 2, with the first 
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image in Figure 2 being a dark view of the same wave dark matter density so that more detail in the 
center of the image may be seen. 

All of these computer generated images of interleaved shells in the wave dark matter density should 
be compared to actual photos of the interleaved shells of elliptical galaxies, such as the image of NGC 
474 in Figure 15. Again, there is no way to know at this time if interleaved shells in wave dark matter 
actually contribute to visible interleaved shells in elliptical galaxies. This is a question that deserves 
further study. Our purpose here is to highlight this intriguing possibility. 

Finally, Figure 16 shows what the wave dark matter density would look like from various angles if 
the wave dark matter glowed white. Note that when viewing the wave dark matter density from along 
the red x-axis that the interleaved shells are not visible. The closer that the point of view is to the x-axis, 
the less pronounced the shells are. This is important because most elliptical galaxies, like the ones in 
Figure 7, do not have visible shells. We should also note that we chose = |^o| to maximize the 
magnitude of the interleaved shells so that they would be easily visible in the computer generated images. 
Choosing |^i| < \Ao\ makes the interleaved shells as subtle as one wants, until they disappear entirely 
when Ai = 0. 

2.2 More Analysis 

Now that we have taken advantage of the fact that a "picture is worth a thousand words" we return to 
equations 21 and 22 and Figure 10 to see how these interleaved shells may be understood in analytical 
terms. For example, we may write the wave dark matter density from equation 22 as the sum of a 
spherically symmetric term, an axially symmetric term, and a third term which we will call the interleaved 
shell density. 

I^WDM = IJ-Spherical + /J'Axial + /^'Shells (23) 

where 

t^Spherical 
1^' Axial 
IJ'SheUs 

When t = 0, a = (j) so that 

IJ-Shells = 2Ao-4i cos(^sin6l • r ■ fu}ofi{r)fuji,iir) = 2AqAi ■ x ■ fujo,o{'^)fuji,i{r). (27) 

It follows that the interleaved shell density, which may be positive or negative, is an odd function. The 
oddness of ^shells explains why the shells are interleaved: When j-isheiis is a local maximum in the 
middle of a shell, it is a local minimum at the antipodal point on the other side of the elliptical galaxy. 
The shells themselves are created when fujoflir) and fu>i,i{r) go in and out of phase with one another, as 
shown in Figure 10. When x is positive, a shell is created when the two functions are roughly in phase 
with one another When x is negative, a shell is created when the two functions are roughly out of phase 
with each other. Note that the two function are in phase with each other three times in Figure 10, and 
that there are three major shells on each side of the origin in the computer generated figures that follow. 

2.3 Conclusion 

Hence we see that wave dark matter might predict interleaved shells in the density of dark matter Since 
interleaved shells are observed in the images of elliptical galaxies, this quaUtative similarity is both 



= Alf^,,o{rf (24) 
= Alsm^e-r^ ■ U^iirf (25) 
= 2AoAicosasm9-r- fojfjfi{r)foj^^i{r). (26) 
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very intriguing and represents another good reason to continue to study wave dark matter as a serious 
candidate for dark matter. 

This global interleaved shell structure in the wave dark matter density exemplified in Figures 2,11, 
12, 13, 14, and 16, while invisible itself, could conceivably play a role in the formation of visible in- 
terleaved shells in some elliptical galaxies, perhaps as exemplified in Figures 1 and 15, through gravity, 
friction, dynamical friction, and other processes. This conjecture is based primarily on the striking qual- 
itative similarities in the previously mentioned images, namely that both sets of images have interleaved 
shells. The basic idea that regions of increased dark matter density could gravitationally attract more 
regular visible matter and hence contribute to visible shells is plausible but is also a very subtle question 
to study. Resolving this conjecture with a high degree of confidence may require sophisticated com- 
puter simulations based on expert modelings of many different astrophysical processes, making this a 
challenging, yet very worthwhile, conjecture to study. 
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A Comment for the Matlab Function Shells.m 



We include one more comment about the matlab function shells.m which produced all of the computer 
generated images of the interleaved shells in the wave dark matter density in this paper. The comment 
is that the wave dark matter density in shells.m was expressed in terms of spherical harmonics which, 
while not needed for the discussion in this paper, is needed if one wants to compute the corresponding 
gravitational potential in the same manner as in [7]. 

In order to expand equations 21 and 22 into spherical harmonics, we need to recall that spherical 
harmonics defined on the unit sphere are actually restrictions of homogeneous polynomials of the same 
degree which are harmonic in R^. Here is a short list of homogeneous harmonic polynomials: degree 
zero: 1, degree one: x, y, z, degree two: — y"^, 2xy, "iz^ — r^, 2xz, 2yz, where = + + z^. 
It is easy to check that these are all harmonic in and hence their restrictions to the unit sphere are 
spherical harmonics of the same degree. 

Since cos asinO = cos(0 — 4>q) sin 6* is a rotated degree one spherical harmonic (the restriction of x 
to the unit sphere) and hence a spherical harmonic itself (taking (/)o to be fixed), the last term of equation 
22 is already in the form of a spherical harmonic times a function of r, as we desire. The first term is as 
well, since it is already a function of r and the zeroth degree spherical harmonic is the constant function 
one. Hence, it is only the middle term that we need to put into the desired form. To do this, we note that 

sin^ = r'^ - ='^r^ + ^(r^ - 3^^), 

since z = r cos 9, which successfully expresses this term as the sum of two terms, each of which is a 
function of r times a spherical harmonic. Hence, 

fiDM ~ Uo{r) + [/2(r)(3z2 - r^) + C/i (r)(r cos a sin ^) (28) 

where 

Uo{r) = Alf^,,o{rf+'^Ayf^,,i{rf 

U2{r) = -3^?/.,i(r)2 

Ui{r) = 2^o^i/^o,o(r)/a;i,i(r). (29) 

These expressions for the wave dark matter potential are the ones used in the matlab function shells.m, 

available at the author's "Wave Dark Matter Web Page" at 
http://www.math.duke.edu/~bray/darkmatter/darkmatter.html. 
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